Collective vortex phases in periodic plus random pinning potential 
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We study theoretically the simultaneous effect of a regular and a random pinning potentials on 
the vortex lattice structure at filling factor one. This structure is determined by a competition 
between the square symmetry of regular pinning array, by the intervortex interaction favoring a 
triangular symmetry, and by the randomness trying to depin vortices from their regular positions. 
Both analytical and molecular-dynamics approaches are used. We construct a phase diagram of the 
system in the plane of regular and random pinning strengths, and determine typical vortex lattice 
defects appearing in the system due to the disorder. We find that the total disordering of the vortex 
lattice can occur either in one or in two steps. For instance, in the limit of weak pinning, a square 
lattice of pinned vortices is destroyed in two steps. First, elastic chains of depinned vortices appear 
in the film, but the vortex lattice as a whole remains still pinned by the underlying square array 
of regular pinning sites. These chains are composed into fractal-like structures. In a second step, 
domains of totally depinned vortices are generated and the vortex lattice depins from regular array. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 

Magnetic and current-carrying properties of supercon- 
ducting films with nanoengineered arrays of periodic 
and quasi-periodic pinning sites attract a lot of atten- 
tion both from the experimenta l 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 and 
theoretical^^'^ii 5 -'^'^ 7 -'^ points of view. The main rea- 
son of this interest is that the regularity in positions of 
pinning sites produces a large increase of the critical cur- 
rent at certain values of the magnetic field which corre- 
spond to the matching between the number of vortices 
and the number of pinning sites in the system. The high- 
est value of the critical current in films with periodic pins 
was obtained for the same concentration of vortices and 
pinning sites. If the pinning potential is strong enough, 
vortices in a film with a square array of pinning sites also 
form a square lattice instead of the triangular one, the 
latter being energetically favorable in the absence of pin- 
ning. This regime corresponds to single-vortex pinning 
when vortices are pinned individually and the flux-line 
lattice behaves not as a collective media. However, if the 
pinning strength becomes of the order of the elastic en- 
ergy of the flux-line lattice, the vortex-vortex interaction 
becomes important, and the triangular symmetry of the 
vortex array is recovered. Thus, at the filling factor 1, 
low values of the pinning strengths favor a deformed tri- 
angular vortex lattice, which is basically depinned from 
the pinning array. At high values of the pinning strength, 
the regular square array of pinned vortices is the lowest 
energy state. There is also an intermediate regime, where 
the vortex lattice is regular, but half of the vortex rows is 
depinned and form a lattice with a symmetry close to the 
triangular one. These three vortex phases are depicted 
in Fig. 1. 

Structural phase transitions due to the tuning of 
the regular pinning strengths were studied theoretically 



within the framework of the London theo^ r 19 : 20 : 21 where 
vortices were treated as point-like objects, i.e., near the 
first critical field. The same regime was considered 
in Refsi 22 ' 23 by using molecular dynamics simulations. 
Refi^i analyzes such transitions near the second criti- 
cal field using the linearized Ginzburg-Landau theory, 
whereas in Refsi 16 ' 17 ! 25 the whole system of the Ginzburg- 
Landau equations was solved numerically. There are also 
applications of the same ideas for the vortex lattices in 
rotating Bose-Einstein condensates of alkali metal atoms, 
where regular pinning potential can be formed and eas- 
ily tuned, in a broad range, by using lasers. At first, 
structural transitions of vortex lattices in this system 
were predicted theoretically in Refsj 26 ' 27 i 28 using Gross- 
Pitaevskii equation and then observed experimentally 29 . 
We also would like to mention a recent experimental 
work 30 on elastic lattices of millimeter-size charged par- 
ticles in a square array of traps which exhibit similar 
physics as the here considered vortex phases. Another 
example is colloidal particles on periodic substrate a 31 ' 32 . 
Thus, the problem of the competition between the sym- 
metry of the underlying pinning array and the lattice of 
repealing particles is important in various fields of mod- 
ern physics. Note that in the case of superconducting 
films, holes which are usually used to create regular pins 
produce a strong pinning potential, which always domi- 
nate the vortex-vortex interaction. Weak enough pinning 
sites can be fabricated, for example, by imposing a pe- 
riodic modulation on the film surface and were recently 
studied theoretically in Refi^. 

In the present paper, we consider theoretically a more 
complicated situation, when there is an additional fac- 
tor, namely, a random pinning potential. Random pin- 
ning is always found in any real physical system; more- 
over randomness cannot be ruled out even in numerical 
simulations, where its effect should be carefully analyzed. 
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Samples potentially important for applications have to be 
very large compared to the period of the regular pinning 
potential, and therefore disorder and defects in lattices 
of vortices or other interacting objects are unavoidable. 
For this reason, understanding of basic properties of such 
systems is of scientific and hopefully of practical impor- 
tance. 



Thus the structure of the vortex lattice in the system 
under consideration is determined by three factors: (i) 
a square array of regular pinning sites, which tries to 
impose its own symmetry on the vortex lattice, (ii) an 
interaction between vortices, which favors a triangular 
configuration of the vortex lattice, (iii) a random pinning 
potential, which attempts to destroy the regularity in the 
vortex positions and to depin them from the square lat- 
tice. We analyze here analytically various kinds of vortex 
lattice defects for a two-dimensional system and deter- 
mine typical values of the random pinning strengths lead- 
ing to a spontaneous generation of these defects. This 
allows us to construct an approximate phase diagram of 
the system in the plane of regular pinning strength and 
random pinning strength. We also perform molecular- 
dynamic simulations in order to model numerically the 
process of disordering of the square vortex array. The 
obtained vortex patterns and the typical values of ran- 
dom pinning strength producing disordering are in good 
agreement with our expectations obtained from the ana- 
lytical treatment. 



Note that the problem studied can be also linked to 
other situations when there is a competition of random- 
ness and regularity, but where the dimensionalities of 
the pinning potential and an elastic media are possi- 
bly different from two-dimensional. These are flux-line 
lattices in layered superconductors, where layers act as 
one-dimensional pinning centers giving rise to a so-called 
intrinsic pinning. Such pinning centers can also be cre- 
ated by twin boundaries. Since there is always disorder 
in real superconductors, one can also have a competition 
between randomness and regularit y 33 ' 34 ! 35 . Other exam- 
ples are charge and spin density waves. 



The paper is organized as follows. In Section II, we 
present the basic formulation of our model and discuss 
our approach of molecular dynamic simulations. In Sec- 
tion III, various kinds of kinks or defects in vortex lat- 
tices are found, their structures are discussed, and their 
energies and sizes are estimated. In Section IV, we char- 
acterize a random potential and explain how to account 
its effect on vortex lattices. Our main results are given 
in Section V, where we construct an approximate phase 
diagram of the system in the plane of regular and random 
pinning strengths. Finally, we conclude in Section VI. 



II. MODEL 



A. Basic formalism 



In our model, we treat vortices as point-like objects 
interacting via a pair-wise potential, which are valid as- 
sumptions as long as the applied magnetic field is much 
lower than the second critical field H C 2 , and the Ginzburg 
Landau parameter is large, k 1. Under such assump- 
tions, the interaction energy of two vortices positioned at 
ri and r 2 can be well approximated by the well-known 
London expression 



H int (ri,r 2 ) 
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(1) 



where Kq is a modified Bessel function. Here and below 
the following dimcnsionlcss variables are used: distances 
are measured in terms of the London penetration depth 
A and energy is measured in units of H 2 A 3 /8ir with H c 
being the thermodynamic critical field. An important 
quantity, which we are going to use, is the interaction 
energy of a regular vortex row with a given vortex situ- 
ated outside of this row. If we put a center of coordinates 
at one of the vortices in the row, and the y-axis is along 
the row, then the interaction energy of the vortex located 
at [x, y) is given by 
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where d is the intervortex distance. Using Fourier trans- 
formation for a modified Bessel function Ko(r) and per- 
forming a summation in the reciprocal space, Eq. (2) can 
be rewritten as 
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Vortices are located in a two-dimensional film with a 
periodic square array of pinning sites. To describe the 
pinning potential of a single site, we use a parabolic func- 
tion: 



U(r) 



-U 1 




(4) 



inside the well, r < a, and the potential is zero outside 
the well, r > a. In this paper, we consider the situation 
when the well radius a is much smaller than the intersite 
distance a, i.e., a <C a. 

In addition to the regular pinning potential there is 
also a random pinning potential in the system; we denote 
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FIG. 1: (Color online) Vortex lattice structure in a super- 
conducting film with weak periodic pinning at filling factor 1. 
Vortices are shown by red (dark gray) filled circles and pin- 
ning sites by black open cirles. Fig. 1(a) corresponds to the 
case of very low pinning with the deformed triangular vortex 
lattice. Fig. 1(b) shows a half-pinned regular phase realized 
at intermediate values of pinning strength. Fig. 1(c) repre- 
sents a pinned vortex lattice, which is energetically favorable 
at higher values of pinning strength. 

its energy by e(r). The total energy of the system is thus 
given by 

E = Hint (r< , tj) + V(&n e ( r *), ( 5 ) 

i,j i,n i 

where and R„ stand for the positions of vortices and 
regular pinning sites, respectively. We also assume that 
the number of vortices is exactly equal to the number of 
pinning sites, i.e., filling factor is one. 

Let us recall first the situation when there is no ran- 
dom potential in the system. In this case, as was ex- 
plained in the Introduction, strong regular pinning favors 
a square vortex lattice symmetry, see Fig. 1(c). At low 
Uq the repulsion between vortices dominates resulting in 
a slightly deformed triangular lattice that has the low- 
est energy, as illustrated in Fig. 1(a). Finally, in the 
intermediate range of Uq, a half-pinned lattice becomes 
the ground state, see Fig. 1(b). In the half-pinned phase, 
an effective pinning potential for depinned vortices is cre- 
ated by their interaction with neighboring rows of pinned 
vortices. It is easy to see that the following expressions 
define the boundaries between the triangular lattice and 
the half-pinned phase, and half-pinned phase and square 
lattice, respectivel y 19 ' 21 

U ~ 2 {E hp ~ Etr) (6) 

U = 2 (E sq - E hp ) , (7) 

where E tr , E} lp , and E sq are the energies of triangular, 
half-pinned and square vortex lattices, without any pin- 
ning, taken per one vortex. In other words, these vortex- 
vortex interaction energies depend only on vortex lattice 
symmetry. We focus on the situation when the strength 
of a regular pinning is not too high, i.e., is close to the 
boundaries defined by Eqs. (6) and (7). In this case the 
vortex lattice can be considered as a collective object and 
we are far from the regime of single- vortex pinning. 



B. Molecular-dynamics simulations 

In addition to an analytical approach, we supplement 
our study with computer experiments. In this Subsec- 
tion, we describe briefly the basic ingredients of the used 
computer simulations. We use molecular-dynamics sim- 
ulations of vortices moving under the action of the forces 
due to the vortex-vortex interaction and the interac- 
tion of vortices with regular and random pinning sites. 
To find the lowest-energy vortex configurations, we per- 
form simulated annealing simulations by numerically in- 
tegrating the overdamped equations of motion (see, e.g., 
Refs, 12 ! 13 - 14 ! 1 ^ 18 ): 

ryv, = U = FT + fT + ff. (8) 

Here, fi is the total force per unit length acting on vortex 
i, f™ and f" p are the forces due to the vortex- vortex 
and vortex-pin interactions, respectively, and f t T is the 
thermal stochastic force. In Eq. (8), r\ is the viscosity, 
which is set here to unity. The force due to the interaction 
of the i-th vortex with other vortices is 

r = -E m^-rji)^ , (9) 

where N v is the number of vortices, K\ is the modified 
Bessel function, and fy = (r, — Tj)/ | — Tj |. The 
modified Bessel function Ki(r) decays exponentially for 
r larger than 1 (A in dimensional units), thus it is safe 
to cut off the (negligible) force for distances larger than 
5. In our calculations, the logarithmic divergence of the 
vortex-vortex interaction forces for r — > is eliminated 
by using a cutoff at distances smaller than 0.1. 

The force due to the interaction of the ith vortex with 
the parabolic pinning potentials (Eq. Q) is: 

N 

^ = ^^)|r i -rW|e(r p -|r,-p« |) r&\ 

(10) 

where N p is the number of pinning sites, f p is the maxi- 
mal pinning force of each potential well, a is the spatial 
range of the pinning potential (for random potentials, 
we define f pr < f p and a r < a as the maximum pinning 
force and radius, correspondingly), O is the Heaviside 
step function, and = (r. t - r[ p) )/ | r l - | . 

The temperature contribution to Eq. ([5]) is modeled by 
a stochastic term obeying the following conditions: 

(flit)) - 0, (11) 

and 

{ff{t)fJ{t'))=2r 1 k B T8 l] 8{t-t'). (12) 

The ground state of a system of vortices is obtained as 
follows. First we set a high temperature to let vortices 
move randomly, then temperature is gradually decreased 
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down to T = 0, thus simulating field-cooled experiments 
(see, e.g. Refs,2L2Z). Note that unlike in RefsJ^ili^ 
we do not consider vortex dynamics driven by an external 
Lorentz force, but we only calculate the ground-energy 
vortex configuration (see, e.g. ) 38 i 39 ). 

In Sec. Ill we discuss the effect of a random poten- 
tial on the square and the half-pinned phases and define 
typical vortex lattice defects induced by the randomness. 



III. KINKS 

Both square and half-pinned phases are periodic in the 
absence of random pinning sites. Translational motion of 
the vortex lattice as a whole by a period of a regular pin- 
ning lattice a does not change the energy of the system. 
Therefore, one can expect that there exist such solutions 
when, in one region of the system, the vortex lattice 
is shifted by a with respect to that in another region. 
There should also be a kink with intermediate value(s) 
for deviations of vortex positions separating these do- 
mains. In this sense, our system resembles a traditional 
sine-Gordon system, which describes the behavior of elas- 
tic media in a periodic potential^ -. For sine-Gordon sys- 
tems, a typical kink has a smooth structure and its length 
depends both on the elastic properties of the media and 
the strength of the periodic potential, and it can be arbi- 
trary large. Kinks are known to play an important role 
in the process of disordering of media 40 . Thus one may 
say that disordering is occurring through a proliferation 
of kinks. Therefore, we will pay special attention to the 
kink structure. 



A. Continuous kinks versus sharp defects 

One can naively assume that there also exist smooth 
kinks for square pinned and half-pinned lattices. How- 
ever, despite the similarity with sine-Gordon systems, our 
system has one specific feature, namely, square and half- 
pinned vortex configurations are not the lowest energy 
states in the absence of pinning. Moreover, they are lo- 
cally unstable, i.e., even infinitesimally small displace- 
ments of vortices lead to the collapse of the lattice into 
a triangular array. In other words, the elasticity the- 
ory does not work properly. This is shown explicitly in 
the Appendix, where we also determine a structure of an 
infinite-lenght kink. We found that it has a "discrete" 
(i.e., not smooth) structure, which is shown in Fig. 2(a). 
This is a fundamental difference between the system at 
hand and conventional sine-Gordon systems. 

Note that in the infinite system this type of defects has 
a divergent energy because it is of infinitely long. There- 
fore, fluctuations generate reconnected kinks of finite sizes 
rather than infinitely long defects, as shown in Fig. 2 (a). 
This point will be explained in detail in Section V. Also 
notice that the fact of appearance of reconnected defects 
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FIG. 2: (Color online) Schematic view of typical defects in a 
regular vortex lattice. Vortices/pinning sites are shown by red 
(dark gray) filled circles/black open cirles. Fig. 2(a) shows 
elastic strings of unpinned vortices in the square pinned lat- 
tice: straight (on the left) and reconnected (on the right). 
Fig. 2(b) presents a row of pinned vortices in the half-pinned 
phase. Fig. 2(c) corresponds to the quasi-one-dimensional 
defect of finite length inside the pinned square lattice. 



is well-known in the condensed matter physics, see e.g. 
Ref. 40 for sine- Gordon and related systems. 

It is important to notice that the string of unpinned 
vortices behaves as an elastic object, vortices inside the 
string are correlated and behave collectively due to their 
interaction. It is easy to see from Fig. 2(a) that such 
a defect can be considered as a nucleus of a half-pinned 
vortex phase. The defect energy density can be there- 
fore estimated as a difference in the energy density be- 
tween the half-pinned and pinned vortex phases. They 
are given, respectively, by 



E. 
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hp 
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(13) 



(14) 



The first term in the numerator in the right-hand side 
(RHS) of Eq. (13) (Eq. (14)) is the energy of the vortex- 
vortex interaction, taken per one vortex, while the second 
term denotes the pinning energy per one vortex. The 
latter terms are different for the half-pinned (Eq. (13)) 
and for pinned (Eq. (14)) phases, since in the second 
case all the vortices are pinned, while in the first case - 
only half. Denominators in the RHS of Eqs. (13) and 
(14) give the area per one vortex. Finally, for the defect 
energy density we obtain 



E 



(i) 



Uq/2 — (E sq — Ehp) 



(15) 



It is important to note that the kink energy vanishes at 
the transition between the square and half-pinned lat- 
tices, i.e., at Uq = 2 (E sq — Ehp) since energies of these 
two phases become equal. 

It is easy to realize that the pinned phase can also be 
disturbed by a domain composed by a region of a de- 
formed triangular vortex lattice. Again such a domain 
is an elastic object, vortices within the domain are de- 
pinned collectively. The region of depinned vortices has 
higher energy density compared to the previously de- 
scribed string-like defect, since the energy of a deformed 
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triangular lattice has a higher energy than the half pinned 
lattice in the relevant parameter region. This energy den- 
sity can be estimated in a similar manner, as the one 
which led us to Eq. (15). For that we need to know the 
energy density of the deformed triangular lattice, which 
is given by»2i 



Etr 
2 ' 



(16) 



where E tr is the vortex-vortex interaction energy in the 
ideal triangular lattice taken per one vortex. In reality, 
this lattice is deformed and it is not ideal. However, as 
it was shown in Refill, the energy of elastic distortions 
is compensated by the pinning energy in the linear ap- 
proximation with respect to vortices displacement, and 
Eq. (16) is correct. Finally, for the density energy of the 
defect we obtain 



E 



(2) 



E, 



(17) 



Now we turn to the description of the half-pinned 
phase. By similarity with the pinned phase, we can see 
that a typical kink in the half pinned phase consists of a 
row of pinned phase. One can also have domains of trian- 
gular phases. The first type of defect is shown schemat- 
ically in Fig. 2(b). The energy densities for these kinks 
can be calculated from the differences in energy densities 
of the nucleating phase and the initial one, as for Eqs. 
(15) and (17). They are given, respectively, by 



E 



(3) 



—Uq/2 + (E sq — Eh p ) 



E 



(4) 



Uq/2 — (E hp — E tr ) 



(18) 



(19) 



which of the kinks has a lower-energy density depends 
on the position in the phase diagram. As again can be 
expected, defects of the first type, shown in Fig. 2(b), 
should be reconnected in finite systems. 



B. Quasi-one-dimensional defects 

We note that a one more type of the defects appears 
in the system for sufficiently larger values of Uq, when 
the potential well for a pinned vortex row is much deeper 
than the effective potential well created for a depinned 
vortex row by surrounding rows of pinned vortices. In 
this limit, the defect considered is a finite row of un- 
pinned vortices, but not reconnected. Such a structure 
can exist due to the fact that in one part of the defect, 
the string of vortices is stretched so that there is one miss- 
ing vortex and the last vortex in this part of the string 
deviates from its equilibrium position by one period of 
the regular lattice a. In another part of the defect, the 
string of vortices is compressed so that there is an ex- 
cess vortex. This defect is shown schematically in Fig. 



2(c) for the most simple case, when it is straight. It 
consists of a finite elastic string of vortices in the same 
row, which are displaced along the direction of the row. 
Thus this object as a whole can be considered as a quasi- 
one- dimensional kink-antikink structure. In Section V, 
by using molecular-dynamic simulations and rather gen- 
eral arguments, we argue that both reconnected chains of 
depinned vortices and quasi-one-dimensional defects are 
limiting cases of wider class of collective fractal-like de- 
fects. In other words, defects considered in this Section 
are not necessarily straight, but straight defects can be 
more easily studied analytically. 

We now calculate typical energy and length of this 
quasi-one-dimensional defect. We assume that positions 
of all other vortices in the system remain unchanged. The 
potential energy of a vortex in the string is given by the 
interaction energy with the fixed vortices, with pinning 
sites, and also with each other. The first contribution to 
the string energy can be found from Eq. (3), where one 
can neglect all the terms except of the ones with m = 0, 
±1 provided that a ~ 1, which results in 



E 



(ID) 



4tt tt ^ ( exp (-a) a . . 2ny k 

— 77 — / -; —, r + - cxp (-27r)cos 

K a \ 1 — exp (—a) tt a 

(20) 

where the summation is performed over positions of all 
vortices in the string y^. Note that from this equation 
one can see that the intervortex interaction favors a shift 
of a/2 of the row in the y direction with respect to the 
lattice making a symmetry closer to the triangular one. 
Note that the potential energy (20) for each vortex has 
a typical sine-Gordon structure. 

The interaction energy for the vortices with pinning 
sites in the row reads 



E. 



(ID) 



E En 



an 



Vh), 



(21) 



where V(y) is defined by Eq. (4) and the index n stands 
for the positions of the pinning sites. For illustrative pur- 
poses, we plot the sum of e[ 1D ^ and E^ D ^ in Fig. 3(a) 
for vortices in the string that are positioned periodically, 
Vk — V + ka. The resulting potential is, of course, a pe- 
riodic function in y. The energy minimum for the string 
is attained for a square configuration due to our choice 
of parameters, i.e., for strong enough regular pinning. 
To show that this potential is of a quasi-one-dimensional 
form, we also plot it in three dimensional form in Fig. 

3(b). 

There is an additional contribution to the energy of the 
vortex, which is due to its interaction with other vortices 
in the string, i.e., the elasticity term. This contribution 
depends on the mutual positions of the vortices. We will 
take into account only the interaction with nearest neigh- 
bors, as is usually done in elasticity theory. Taking into 
account Eq. (2) and using a Taylor expansion for the 
modified Bessel function, we obtain the following expres- 
sion for the energy of the whole string of vortices 
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There are Im/a vortices along this defect and they are 
dcpinned from their lowest energy positions. Each vor- 
tex gains some energy by being depinned. If the width 
of the potential well a is much smaller than the inter- 
site distance a, the gain is given by Uq, as follows from 
Eq. (23). At the same time, each vortex also loses some 
amount of energy due to the fact that it no longer forms a 
square configuration with vortices from the surrounding 
fixed rows, as can be seen from Eq. (20). Thus, the total 
energy increase for the string of Iid/® depinned vortices 
is given approximately by 



(ID) 



E. 



(ID) 



4tt 



Iid/o. U je 



The elastic energy of the string E^ 1D ^ can be easily esti- 
mated from Eq. (23) as 



-2tt 



(25) 
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K (a) 
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(26) 



If we now minimize the total defect's energy E^ 1D ^ with 
respect to its length Im, we obtain 



27Tfl 4 



D 



K {a) + -^\. (27) 



This leads to the following result for the defect energy 



FIG. 3: Potential energy profile for the vortex row in the 
square lattice along one of the principal axis, provided that 
positions of vortices in the row are fixed with respect to 
each other and all other vortices are pinned. Fig 3(a) cor- 
responds to the dependence on one coordinate, whereas Fig. 
3(b) shows a dependence on both coordinates, which is of 
quasi-one-dimensional form. 
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where Auk is the displacement of kth vortex with respect 
to the (fc — l)th one. As usual, we can replace Auk/a 
by du/dy and also switch from the summation to the 
integration in the RHS of Eq. (22)i£: 



E- 
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dy, (23) 



where the integration is performed along the length of 
the defect. The total energy of the defect is thus given 
by the sum of three contributions: 



E^ = E\ 



(ID) 
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(24) 



which can be estimated in a simple manner. We assume 
that a typical kink-antikink structure has a length 1\d- 
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(ID) 
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which has a nontrivial dependence on Uq. As we will 
see below, in the limit of weak pinning, quasi- ID defects 
have higher energies compared to previously analyzed 
ones. This can be explained by noticing that vortices 
in the string are situated neither exactly in a potential 
well induced by pinning sites nor in that created by the 
neighboring rows of vortices. 

Note that a quasi- ID defect can have a slightly dif- 
ferent structure, when the maximum deviation of vor- 
tex in the string from its equilibrium position is given 
by a/2 but not by a. In this case, the position of the 
most strongly displaced vortex in the row is stabilized by 
the interaction with neighboring rows of vortices and not 
by a pinning site. Our estimates for length and energy 
of quasi-one-dimensional defect, given by Eqs. (27) and 
(28), however remain applicable also for such a defect. 



IV. RANDOM POTENTIAL 

In this Section, we characterize the random potential 
and describe how to deal with it in our problem. We 
assume that the random potential is created by randomly 
distributed pinning wells with size cr r < a, a and depth 
U r < Uq. The concentration of random sites is denoted 
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by n r , and we also assume that it is higher or equal to 
the concentration of regular pinning sites, n r > a~ 2 . 

It is well-known that the actual pinning energy due to 
the random pinning is produced by fluctuations of the 
random potential in a given region of the system, i.e. by 
bunches of pins, and not by pinning potential averaged 
value. If we have a region of area So 3> n^ 1 , the distri- 
bution probability P(N) for having N random pinning 
sites within this region can be well described by a normal 
distribution 



P(N) = 



1 



1 



/27r >/n r So 



exp 



(N - n r S Q ) 2 
2n r So 



(29) 



As can be expected, the averaged value of N is given by 



(N)=n r S , 
whereas the mean-square deviation reads 



(N 2 ) - (AO 2 - n r S . 



(30) 



(31) 



The probability to have any number of random pinning 
sites larger than some definite value No within the region 
So can be found by integration, P(N)dN, and this 
integral can be estimated as (provided that No — n r So > 

V n rS ) 



P(N)dN 



1 



V n rSo 



No 



/2tt N q - n r S Q 



exp 



(N - n r Sof 



2n r So 



(32) 

The probability for this event not to occur is, of course, 
1 — f Ng P(N)dN. If we now have a larger region S > So, 
the probability that, within S, one cannot find a re- 
gion of area So with number of random pinning sites 

( \ S I So 
l — f No P(N)dN) 

Finally, the probability to find such a region is 1 — 
1 — P(N)dN^j . From this analysis and from 

Eq. (32) it is easy to see that if S ~ So and \Nq — n r So | ~ 
y/n r So, this probability is very close to 1. This means 
that within the domain of area S slightly larger than 
So we can always find a bunch of pines, such that the 
pinning energy around this bunch (within region of area 
So) deviates by y/n r So from its average value. This il- 
lustrates that the effective pinning energy within some 
region can be closely approximated by mean-square de- 
viation of random pinning potential within this region, 
U r \/n r So- This conclusion is well-known in the physics 
of disordered flux- line lattices 41 . The famous Larkin- 
Ovchinnikov theory is also based on the idea that the 
mean-square deviation of random pinning energy is bal- 
anced with the energy of lattice distortions within some 
domain. 

Each kind of defects described in Section IV has its own 
typical energy and sizes. In order to estimate the ran- 
dom pinning strength leading to a proliferation of such 
defects one can equate the defect energy to the mean- 
square deviation of the random pinning potential within 



the characteristic area of the defect. By doing this, one 
can study the phase diagram of the system in the plane 
of random and regular pinning sites strengths. 



V. PHASE DIAGRAM 

In the absence of random pinning, there are two phase 
boundaries separating the pinned square, half pinned and 
deformed triangular lattices, as given by Eqs. (6) and 
(7). We now consider the effect of the random pinning 
potential. Note that a deformed triangular vortex lat- 
tice is aperiodic even in the absence of random pinning, 
since periods of vortex and regular pinning sites lattices 
are incommensurate. Therefore we consider disordering 
of square pinned and half-pinned vortex lattices, which 
have a regular structure. As was already explained, the 
destruction of the order occurs via a proliferation of de- 
fects. Using the tools developed in Sees. Ill and IV, we 
can estimate the typical random pinning strengths lead- 
ing to a generation of different defects, identified in Sec- 
tion III. We do not expect that this will give us the full 
phase diagram, since we have found defects only in some 
simple and limiting cases. In order to obtain a deeper in- 
sight to the problem we will also use molecular dynamic 
simulations. 

Let's start from the case of a square pinned lattice. The 
lowest energy defect in such a situation, among those we 
have found, is constructed from the elastic string of col- 
lectively depinned vortices, as explained in Section III. 
This row has to be reconnected because otherwise its 
energy is infinitely large 40 (therefore, one should pay a 
special attention to this point both in numerical and real 
experiments on finite arrays, since there it is easy to get 
just straight defects - their energies remain finite because 
the system sizes are finite). The radius of such a defect 
cannot be too small, since the system is discrete. There- 
fore, the lowest energy defect has a radius Rd of the order 
of the intersite distance, Rd ~ a, see Fig. 2 (a). Let us 
write Rd as ar<j, where the dimensionless factor is of 
the order of one, r<j ~ 1. The energy of such a defect can 
be found from Eq. (15): 



E k 



ink 



27rr d [U /2 - (E sq - E h 



)}. (33) 
To find the 



The area of this structure is just ira 2 
random pinning strength leading to a proliferation of this 
particular kind of defects, we equate the energy of the 
defect -Efci n fc to the mean-square deviation of the random 
pinning energy within the area of the defect U r ^/n r 'Ka 1 r 2 d 
and we find 



[Uo - 2 (E sq - E hp )} 



(34) 



From Eq. (34) we see that the U r required to destroy 
the order in the vortex positions in the square pinned 
phase vanishes at Uq = 2 (E sq — E^p), i.e., in the point 
of the transition between the square pinned and half- 
pinned lattices in the absence of disorder. This means 
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FIG. 4: Schematic phase diagram of the vortex lattice in 
the plane of the periodic pinning strength Uo and the random 
pinning strength U r . 

that arbitrary weak disorder destroys the regularity in a 
square vortex lattice in the vicinity of this point, since 
the energies of pinned and half pinned lattices are equal 
to each other there. 

We now start to construct a phase diagram of the sys- 
tem in the J7o - U r plane for the same concentration of 
random and regular pins, i.e., n r a 2 = 1. The bound- 
ary separating pinned and a mixture of pinned and half- 
pinned lattices is shown by line 1 in Fig. 4 at a = 1 
and k = 5. This mixed phase still can be considered as 
pinned, as a whole, by the regular pinning sites array, 
since there is a strong correlation between the positions 
of regular sites and the vortices. 

Let us now switch to the half-pinned phase. As was 
explained in Section III, the lowest energy defect in this 
case corresponds to the row of pinned vortices, which 
recovers locally square array, as in Fig. 2(b). This row 
again has to be reconnected. By repeating the derivation 
which led us to Eq. (34) , we obtain the following criterion 
for U r favoring proliferation of such defects 



U '~\ 2 t-^o + 2 (E sq - E hp )} . (35) 

It corresponds to the curve 2 in the phase diagram which 
is basically a mirror image of the boundary given by Eq. 
(34). It is interesting to note that disorder in this case 
leads to an unexpected situation, namely, it forces vor- 
tices within some regions to be pinned by the regular 
array of pinning sites. Also note that if we follow along 
curves 1 and 2, we see that the minimal random pin- 
ning strength required for the partial disordering of the 
vortex lattice is a non-monotonic function of Uq. This 
is again a nontrivial result, since one could expect that 
the higher Uo, the higher must be U r , which destroys 
the order in the vortex positions. However, this is not 
exactly the case, since there is one more factor, i.e., the 



intervortex interaction, which helps the random pinning 
sites to destroy order. Thus, at a certain value of Uo it 
is energetically favorable for the system to decrease its 
energy not by introducing vortex lattice defects, but by 
making a structural phase transition of the vortex lattice, 
on large scales, into a half-pinned regular phase. 

The half-pinned phase can also be destroyed via elastic 
domains of the deformed triangular lattice; the energy 
density for the corresponding defect is given by Eq. (19). 
Thus, the random energy U r generating such domains 
can be estimated as 



U r ~ J 2 [U - 2 (E hp - E tr )] . (36) 

y Tly Ct 

It is shown as curve 3 in Fig. 4, which separates regions 
of half-pinned and deformed triangular, i.e., completely 
disordered, phase. Let us stress that proliferation of such 
domains leads to a complete depinning of the vortex lat- 
tice from the underlying square array of pinning sites. 
Curves 2 and 3 intersect in some point and they bound 
a region in the phase diagram where the half pinned vor- 
tex lattice is stabilized. Of course, this does not mean 
that there are no defects at all in this region (or in the 
region of the square pinned lattice) , but this implies that 
their concentration is low. One should also understand 
that the state, which we call a mixture of pinned and 
half pinned phases, contains domains with the deformed 
triangular lattice as well, but again the concentration of 
such domains is low. There can be other types of defects, 
like vacancies, when one vortex is redistributed to the po- 
sition which is very far, in terms of lattice period, from its 
initial position, however, the probability of such events 
is low, since we are in the regime of collective pinning. 
Therefore, such events do not contribute significantly to 
the process of vortex lattice disordering. 

Now we discuss how a mixture of pinned and half 
pinned phases is finally completely depinned from the 
regular pinning array. For that we first estimate U r gen- 
erating domains of triangular lattice inside the square 
one: 



U r ~2J ^[U -(E sg -E tr )], (37) 

y TlyCl 

which is depicted as curve 4 in Fig. 4. In reality, before 
such domains appear in the square phase, it is already 
disturbed by rows of the half-pinned lattice, so that curve 
4 does not describe properly such a transition. The same 
is true for the curve 3 which describes the transition from 
the half pinned to the deformed triangular lattice. In 
other words, curves 3 and 4 must intersect. 

From the obtained phase diagram we see that the dis- 
ordering of the square pinned lattice in the limit of weak 
pinning always occurs in two steps. First, some rows of 
vortices depin and form reconnected elastic chains; their 
concentration increases with the increase of random pin- 
ning strength U r - At higher values of U r , elastic domains 
of deformed triangular lattice proliferate in the system 
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leading to a complete disordering of vortex array. If we 
start from the half pinned lattice, two scenarios are pos- 
sible. If the regular pinning strength Uq is relatively low, 
then the order is destroyed in one step due to the pro- 
liferation of domains of the deformed triangular lattice. 
If Uq is relatively high, again a two-step disordering can 
take place. As a first step, an admixture of the square 
pinned lattice appear in the system, and as a second step, 
domains of a deformed triangular lattice depin the vor- 
tex lattice from the regular pinning array. Note that in 
experimental work 30 on macroscopic charged balls on a 
square array of traps, the formation of domains in the 
lattice of these balls was also observed. Although there 
was no artificial randomness in the experimental setup, 
the disorder was certainly unavoidable. Since energies of 
three relevant vortex phases, namely, triangular, pinned, 
and half-pinned ones, are quite close to each other, do- 
mains formation can be observed in such experiments. 

In Section III, we also have described a quasi-lD type 
of defects, which destroy the order in vortex positions 
within the same row, see Fig. 2(c). The typical area of 
such a defect is ~ ?ir>a, where lm is given by Eq. (27). 
The energy of the defect is defined in Eq. (28). Thus, for 
U r leading to the generation of quasi-lD defects we get 



U r 



27T 
: 72 



U - 



4tt 



3/4 r Kl {ay 1/4 



K (a) 



(38) 

Our estimates show that, for realistic values of the pa- 
rameters, U r could become larger than the values cor- 
responding to the phase boundary depicted by lines 2 
and 3. This conclusion is due to the fact that quasi-lD 
defects represent an intermediate case between the col- 
lective pinning of vortices and a regime of single-vortex 
pinning at high pinning strengths. The vortex lattice in 
this case is still an elastic media, with the intervortex 
interaction being important in the total balance of the 
energies. The phase boundary corresponding to such a 
transition is plotted in Fig. 4 as line 5. Since quasi- ID 
defects destroy the order in the vortex positions in the 
same row, a mixture of pinned and half-pinned phases in 
this case is also characterized by the absence of such an 
order in the rows of unpinned vortices; moreover, chains 
of unpinned vortices are now not necessarily reconnected. 
In this sense, the difference between the quasi-lD disor- 
dered phase and a mixture of pinned and half-pinned 
states is smeared out. 

One can expect that the final regime with sufficiently 
high Uq corresponds to zero-dimensional defects, when 
vortices are depinned individually that is known as the 
single-vortex pinning regime. Thus, the dimensionality 
of typical defects leading to the destruction of the order 
is decreasing continuously from 2 to when the regular 
pinning strength Uq is increasing. At low Uq, regularity in 
the vortex positions is destroyed by the appearance of re- 
connected elastic chains of collectively-depinned vortices 
that are ordered within each chain. At higher Uq, these 
chains start to be open, they should become straighter 



and, therefore, one can expect neighboring chains to con- 
nect due to the lack of free space. If we again increase Uq, 
chains become quasi-one dimensional elastic rows of finite 
length, as we have shown schematically in Fig. 2(c). Fi- 
nally, we are in a single-pinning regime, when each vortex 
is pinned individually, that can be interpreted as chain 
shrinking towards a single lattice period, since the chain 
length given by Eq. (27) decreases with the increase of 
regular pinning strength. 

So far only limiting cases of the first stage of vortex lat- 
tice disordering have been considered, namely, when we 
have: (i) reconnected chains, i.e., two-dimensional col- 
lective defects, (ii) quasi-one-dimensional stripe-like col- 
lective defects, (iii) quasi-zero-dimensional individual de- 
fects in the regime of a single-vortex pinning. To have a 
deeper insight into the problem and to understand better 
how these defects evolve, we perform molecular-dynamics 
simulations for the regime intermediate between the weak 
and strong pinning regimes. In Figs. 5-8 we present 
typical vortex patterns corresponding to an increasing 
strength of the random pinning potential U r with a con- 
stant concentration of random sites, which is set equal 
to the concentration of regular pins. The initial state 
of the system (i.e., without random pinning) is a square 
lattice of vortices pinned on regular pins. A simulation 
region with the sizes 20 x 20 A contains an array of 20 x 20 
regular pinning arrays and 400 vortices. We use periodic 
boundary conditions at the boundaries of the simulation 
region. The value of Uq is 0.045 x 2tt/k 2 , whereas the 
intervortex distance is 1 in units of A. This value of Uq 
is several times larger than that required for the transi- 
tion from the square phase to the half-pinned state, i.e., 
we are indeed in the intermediate regime, where both 
reconnected and open chains of depinned vortices are ex- 
pected to appear in the system with increasing U r . This 
is exactly what we see in Fig. 5 at U r = 0.72£/o, which 
contains both types of chains, and still there is no do- 
mains of a deformed triangular lattice. In general, these 
elastic chains are packed into finite fractal-like structures 
with a quasi-self-similar topology and dimensionality be- 
tween 1 and 2. The building blocks of these fractals are 
represented by parts of rectangles of discrete sizes (i.e., 
commensurate with the period of the pinning array). If 
we increase U r , as in Fig. 6 for U r — 0.85£/o, domains of 
deformed triangular lattice start to appear in the system, 
and thus they coexist with chains of unpinned vortices. 
However, these domains somehow suppress the networks 
of depinned vortices by simply cutting them, since there 
is a competition between these two kinds of defects for 
a free space in the system. Further increasing U r leads 
to the increase of the concentration of domains of totally 
depinned vortices, as indicated in Fig. 7 at U r — 1.28f7o- 
It is interesting to note that domains of pinned square 
vortex lattice turn out to be rather robust with respect 
to the random pinning; some of them survive up to quite 
large values of U r , several times greater than that corre- 
sponding to the situation shown in Fig. 7, with the frac- 
tion of pinned vortices decreasing slowly with increasing 
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FIG. 5: (Color online) Vortex pattern in the film with regular 
square and random pinning arrays for U r = 0.72(7o, obtained 
in molecular-dynamics simulations. Dashed lines are guides 
for eyes indicating positions of chains of unpinned vortices. 
Irregular blue (light gray) spots represent positions of random 
pins, red (dark gray) filled circles show positions of vortices, 
and regular black open circles correspond to periodic pins. 
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FIG. 6: (Color online) Vortex pattern in the film with regular 
square and random pinning arrays for U r = 0.85[/o, obtained 
in molecular-dynamics simulations. Dashed lines are guides 
for eyes indicating positions of chains of unpinned vortices. 
Dash-dotted lines show borders of domains of deformed tri- 
angular lattice which is depinned from the underlying regular 
pinning array. Irregular blue (light gray) spots represent po- 
sitions of random pins, red (dark gray) filled circles show po- 
sitions of vortices, and regular black open circles correspond 
to periodic pins. 




U r . Fig. 8 shows the final state with quite high value of 
U r = 5.3C/o, when almost all the vortices are depinned 
from the regular pinning array and there is no correlation 
between the positions of the regular pins and vortices, but 
there is a pronounced correlation between the positions 
of the vortices and random pins. Note that typical val- 
ues of U r leading to the generation of defects are in good 
qualitative agreement with our analytical predictions. 

We also performed simulations for even higher Uq. We 
found that indeed the initial defects turn out to be quasi- 
one-dimensional stripes with the characteristic length of 
few lattice periods, rather than long-length and branched 
structures. Other defects found correspond to higher val- 
ues of U r , and they are more similar to a single-particle 
imperfections than to regions of elastically distorted me- 
dia, which are constructed of many interacting particles. 
In fact, these defects correspond to domains of totally 
depinned vortices described above, in the limiting case 
of a single-vortex pinning. Finally, single-vortex defects 
and stripe-like defects also match, when the length of the 
stripes tends to a one lattice period. A general qualita- 
tive picture of the evolution of defects is presented in Fig. 
9. 

Thus, the molecular-dynamics simulations enable us 
to relate defects of various dimensionalities and to better 
understand their evolution. 



VI. CONCLUSIONS 

In this paper, we studied theoretically the effect of dis- 
order on vortex lattices in the presence of a regular square 
array of pinning sites. We considered a two-dimensional 
system with the same concentration of regular pinning 
sites and vortices, i.e., with filling factor 1. We were 
mostly focused on a weak pinning regime, when vortex 
lattice demonstrates collective properties due to the mu- 
tual repealing of vortices. The resulting structure of the 
vortex lattice is determined by the competition between 
the square symmetry of the regular pinning array, by 
the triangular symmetry imposed by the vortex-vortex 
interaction, and by the action of randomness trying to 
destroy the order in vortex positions. Typical defects of 
vortex lattices have been identified for some limiting and 
therefore simple cases and their energies and sizes have 
been estimated analytically. Molecular-dynamic simula- 
tions were also used for exploring regions on phase di- 
agram, which are difficult to approach in an analytical 
way. Topological defects found are of various dimension- 
alities including fractional ones having fractal-like struc- 
ture. We were able to predict analytically two limiting 
cases of these fractal-like structures: a reconnected chain 
of unpinned vortices and a straight stripe-like defects, 
while molecular dynamic simulations allowed us to re- 
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FIG. 7: (Color online) Vortex pattern in the film with regular 
square and random pinning arrays for U r = 1.28(7o, obtained 
in molecular-dynamics simulations. Dashed lines are guides 
for eyes indicating positions of chains of unpinned vortices. 
Dash-dotted lines show borders of domains of deformed tri- 
angular lattice which is depinned from the underlying regular 
pinning array. Irregular blue (light gray) spots represent po- 
sitions of random pins, red (dark gray) filled circles show po- 
sitions of vortices, and regular black open circles correspond 
to periodic pins. 



veal the connection between these limits. 

We then constructed an approximate phase diagram of 
the system in the plane of random and regular pinning 
strengths. We found that disordering of the regularity 
in vortex positions can occur either in one or two steps 
depending on the strength of the regular pinning poten- 
tial. In the former case, domains of deformed triangular 
lattice start to appear spontaneously in the film. In the 
latter case, in the beginning, reconnected chains of un- 
pinned vortices are generated inside a square lattice (or 
the opposite: chains of pinned vortices appear inside the 
half- pinned phase) , and only after that again domains of 
deformed triangular lattice depin the vortex lattice as a 
whole from the regular pinning array. Both chains and 
domains are elastic defects with interacting vortices de- 
pinned collectively. 

We also analyze an intermediate regime for values of 
regular pinning strength, which are between the regimes 
of collective and single-vortex pinning. In this case, 
the first step of disordering occurs by proliferation of 
opened, not reconnected, elastic chains of depinned vor- 
tices, which can form long fractal-like defects with com- 
plicated quasi-self-similar topology. Fractals appear due 
to the process of opening of reconnected chains and their 
straightening, when different chains start to connect to 
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FIG. 8: (Color online) Vortex pattern in the film with regular 
square and random pinning arrays for U r = 5.3Uo, obtained 
in molecular-dynamics simulations. Irregular blue (light gray) 
spots represent positions of random pins, red (dark gray) filled 
circles show positions of vortices, and regular black open cir- 
cles correspond to periodic pins. 
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FIG. 9: General picture of the evolution of disorder-induced 
defects for square pinned vortex lattice. Lower branch cor- 
responds to the elastic strings of collectively depinned vor- 
tices. First, they appear as reconnected chains. Then chains 
break and straigthen. As a result, they start to overlap and 
form fractal-like clusters. Finally, they decay to quasi-one- 
dimensional stripes, which shrink gradually to quasi-zero- 
dimensional individual defects corresponding to the single- 
pinning regime. Upper branch shows the evolution of higher- 
energy elastic defects, which form domains of vortices collec- 
tively depinned from regular pins. These domains break into 
uncorrelated blocks of smaller sizes and finally also shrink to 
quasi-zero dimensional defects thus matching with the lower 
branch of defects. 
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each other. At the same time, the average length of these 
defects is going to decrease with the increase of regular 
pinning strength due to a periodic potential sharpening, 
so that there is a competition between two tendencies. 
As a result, chains become straighter and relatively short 
thus forming quasi-one-dimensional stripes. Finally, for 
high values of periodic pinning strength, they shrink con- 
tinuously towards a quasi-zero-dimensional defects corre- 
sponding to a single-pinning regime. The latter regime is 
characterized by one-step disordering, i.e., the width of 
the region on the phase diagram with coexistent chains of 
unpinned vortices and domains of totally depinned vor- 
tices has to vanish with the growth of regular pinning 
strength, since these two types of defects also match. 
Thus, we can understand on the same footing all the 
regimes of disordering starting from a collective pinning 
to a single- vortex pinning. 

Our results are applicable not only to vortices in su- 
perconductors, but also for a broader class of physical 
systems containing a large number of repealing objects 



and underlying system of traps: for instance, vortices in 
ultracold gases in optical lattices, colloids, macroscopic 
elastic balls interacting via Coulomb forces, charge and 
spin density waves. 
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VIII. APPENDIX 

Let us assume that for a square pinned lattice, an 
infinite kink can have a smooth structure and is di- 
rected along the y axis {x and y axes are directed like in 
Fig. 1(c)), thus separating two regions of vortex lattices 
shifted by a in y direction with respect to each other, as 
for usual sine-Gordon systems^. We also denote kink's 
length along x as I2D {I2D ^ a). There are ho/a ^> 1 
vortex rows within the kink's length and we denote devi- 
ations of vortices in the n's row along the y direction as 
u n . We also restrict ourselves to lattices with intervortex 
distances of the order of the penetration depth A (a ~ 1 
in dimensionless units). In this case, the dominant con- 
tribution to the interaction energy of a given vortex with 
all other rows of vortices is provided by harmonics with 
m = 0, ±1, as can be seen from the exponent in the RHS 
of Eq. (3). This energy is then given by 



(27T) 
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Now we can expand the RHS of Eq. (Al) in terms of 
Au n /a. To find the energy of the kink, we have to sum 
the contributions from all the rows within the kink. But 
instead it is convenient to switch from summation to inte- 
gration with the simultaneous exchange of Au n /a by 4^, 
i.e., to use a continuous limit and introduce smoothly 
varying deformation field of the vortex lattice. This 
method is applicable for sine-Gordon systems^, and we 
repeat the same arguments for the system at hand. Fi- 
nally, for the u dependent part of the vortex-vortex en- 
ergy of the kink (per length of the kink in y direction), 
we have: 
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The derivative 4^ can be estimated as 4^ « -p 2 - , and this 

ax ax L2D 



leads us to the following expression 
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There is also a contribution to the kink energy coming 
from the regular pinning potential. It appears due to the 
fact that vortex rows within this kink are depinned. In 
the limit of small potential wells, a -C a, this energy (per 
length of the kink in y direction) is given by 



TP 



(vp) 
kink 



2D 
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For the total energy of the kink we thus have 



TP - TP^ XPW ~ ( 2t „-2tt 1 ,rj l 2D 

f% <'2D a 

If we now try to minimize the kink energy given by 
Eq. (A5) with respect to kink's length, we immediately 
see that I2D = is the lowest energy solution. Of course, 
in reality I2D cannot be shorter than a, and I2D = a is 
a true solution. This conclusion is based on the physical 
fact that a square lattice of free interacting vortices is 
unstable in the absence of pinning, leading to a negative- 
ness of " elastic constant" , in contrast to the sine-Gordon 
and related systems^. 



